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تحلیل قابلیت اعتماد سازه‌ها براساس روش تر کیبی کوچک‌ترین مربعات دستگاه بردارهای پشتیبان و شبیه‌سازی 
مونت کارلوی پیشرفته * 


(یادداشت پژوهشی) 


0( 0 
محمدمهدی مجاهد بهروز احمدی ندوشن 


چکیده بانوجه به کوچک بودن احتمال حرابی سازه‌هاء تحلیل قابلیت اعتماد سازه‌ها هزین زمانی بالایی را به‌همراه نحواهد داشست. در این 
مقاله به‌منظور کاهش زمان محاسبات, الگوریتمی براساس ت رکیب روش رگراسیونی کوچک‌ترین مربعات دستگاه بردارهای پشتیبان و دو 
روش شبیه‌سازی موئت کارلوی پیشرفته: نمونه‌برداری بااهمیت و ابرمکعب لانینی, ارائه شده است. با اراهُ دو مثال قاب و یک مثال حرپا 
کارایی الکوریتم‌های پیشنهادی مورد ارزیابی قرا رکرفته است. تنایج حاصل نشان می‌دهد که روش پیشنهادی می‌تواند احتمال خحرابی را 


به‌نحویی تحمین زند و زمال محاسبات در مقایسه با دیکر روش‌های ارائه‌شده در سال‌های احیر بسیا ر کمتر می‌باشد. 


واژه‌های کلیدی قابلیت اعتماد سازه‌ها؛ روش‌های پیشرفتةٌ شبیه‌سازی مونت‌کارلو دستگاه بردارهای پشتیبان» احتمال خرابی. 
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" تاریخ دریافت مقاله ٩۱/۱۲/۳‏ و تاریخ پذیرش آن ٩۲/۱۱/۹‏ می‌باشد. 


(۱) نویسنده مسئول, دانشجوی کارشناسی ارشد سازه دانشکده عمران. پردیس فنی و مهندسی. دانشگاه یزد. 
حمع. هصرع )م2 رمک وه :لنمم۴ 
(۲) دانشیار دانشکده عمران» پردیس فنی و مهندسی. دانشگاه یزد. 


مقد مه 

قابلیت اعتماد. موضوعی بااهمیت در تحلیل و طراحی 
سازه‌ها می‌باکتة.بیخت ای تحلیار فابلیعت اعتمتاه 
سازه‌ها. محاسبهٌ تخمین احتمال خرابی ناشی از عدم 
قطعیت پارامترها خواهد بود. به‌دلیل ماهیت غیرقطعی 
پارامترهای سازه و در دسترس نبودن تابع حالت حدی 
به فرم ریاضی» محاسبة ایمنی مطلق غیرممکن می‌باشد 
[1]. 

باتوجه به مشکلات روش‌های تحلیلی و 
شبیه‌سازی. از اوایل ده ۸۰ میلادی روش‌های سطح 
پاسخ به‌منظور کاهش حجم محاسبات و درنتیجه کاهش 
هزین زمانی ناشی از روش‌های المان محدود به‌کار 
گرفته شدند. ایدهُ اساسی این روش‌هاء. جایگزینی توابع 
حالت حدی با توابع سطح پاسخ ساده‌تر می‌باشد. با 
تخمین تابع عملکرد مناسب. می‌توان تحلیل قابلیت 
اعتماد سازه‌ها را به‌درستی انجام داد [2-4]. 

در دو دهه اخیر با پیشرفت تحقیقات در زمينهة 
هوش مصنوعیء روش شبکةُ عصبی مصنوعی به‌منظور 
تقریب تابع عملکرد به‌سرعت جایگزین روش‌های 
پیشین شد. کاربرد و جوانب مختلف روش شبکه 
عصبی در بسیاری از مقالات مطرح شده است [5-7]. 
شبکهة عصبی با تمام مزیت‌هایی که نسبت به روش‌های 
کلاسیک دارد. در بسیاری از مسائل سازه‌ای با تعداد 
متغیرهای زیاد می‌تواند هزينة زمانی نسبتاً زیادی» برای 
محاسبات دربر داشته باشد [8]. 

در همین راستا برای رفع مشکلات شبکة عصبی, 
روش رگراسیونی دستگاه بردار پشتیبان ( 5020۲ 
(۱۷]20:06)5۷ ۷۵۵۲۵۶۲) به‌عنوان یک روش یاد گیری 
هوشمند برای به‌دست آوردن مدل رگراسیونی مناسب؛ 
به‌کار گرفته شده است. واپنیک برای اولین بار در سال 
و توا فا سای ور سای سیواة 
داد [9]. روش 5٩۷‏ برای تقریب توابع عملکرد: با 
افزایش تعداد داده‌های آموزشی. کارایی مناسبی از حود 
نشان نمی‌دهد لذا کارایی آن برای سازه‌هایی که دارای 
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تحلیل قابلیت اعتماد سازه‌ها براساس روش ترکیب یک وچک‌ترین.. 


متغیرهای زیادی هستند. کمتر می‌شود و هزينة زمانی 
زیادی را برای آموزش مدل صرف می‌کند. برای رفع 
مشکلات ناشی از این روش, بای گانچن در سال ۲۰۰۹ 
به بررسی کاربردهای روش رگراسیونی کوجک‌ترین 
مربعات دستگاه بردارهای پشتیبان ( ۹002765 1.6251 
185۷ مصنطم۱2 ۷۵۵۲۵۲ 9۷۵00۲۲) و ترکیب 1 
با روش نمونه‌برداری مونت‌کارلوی معمولی پرداخت و 
با ارائ مثال‌های عددی نشان داد که این روش درمقایسه 
با روش ٩۷۷۲‏ از کارایی بالاتری از لحاظ کاهش هزينهة 
زمانی برخوردار است [10]. 

باتوجه به‌اینکه کاربرد روش 1,۹8۷ اخیرا 
به‌عنوان نسخهٌ بهبودیافته‌ای از ٩۷۷‏ در سایر رشته‌های 
علوم و مهندسی بسیار مورد توجه قرار گرفته است 
[10-13] در تعیین قابلیت اعتماد سازه‌ها بسیار محدود 
است و هم‌چنین باتوجه به مشکلات روش‌های نوین 
دیگر برای کاهش زمان محاسبات؛ لذا در این مقاله از 
ارزو انشفاهه کته اتب زما تاک اناد تعیل 
براساس روش شبیه‌سازی مونتکارلوی معمولی زیاد 
شم باه لا جر اس مقالته پاسفاهه از وون‌هبای 
نمونه‌برداری بااهمیت و ابرمکعب لاتینی و براساس 
اطلاعات اضافی مسئله. نواحی انتخابی نمونه‌ها در 
قسمت‌های مهم و موثر فضای متغیرها محدود شده و 
بدین ترتیب تعداد تحلیل‌های مورد نیاز و درنتیجه زمان 
کل محاسبات نسبت به روش شبیه‌سازی مونت‌کارلوی 
معمولی کاهش داده می‌شود. 

درادامه باتوجه به اصولی که هر کدام از روش‌های 
پیشنهادی در کاهش زمان محاسبات دارند. به بررسی 
روش‌های نمونه‌برداری بااهمیت و ابرمکعب لاتینی 
درجهت تولید نمونه‌ها برای تحلیل سازه‌ها و سپس به 
فرمول‌بندی 1,5۷1 پرداخته می‌شود. درنهایت با ارائه 
مها تاو شا ز‌های اسکتی رد فا و یک ریت ): 
کارایی بالای الگوریتم‌های پیشنهادی درمقایسه با 
سایر روش‌های نوین ارائه‌شده در سالیان اخیر جهت 
تحلیل قابلیت اعتماد این سازه‌ها نشان داده شده 


شاج بیست و نهم. شماره یک» ۱۳۹۹ 


محمد مهدی مجاهد - بهرو ز احمدی دوشن 
تاه 


روش نمونه‌برداری بااهمیت 
((157۷۲) ۱۷۲۵۵۵0 عصتامر‌صجه م‌صم) همه 
روش نمونه‌گیری بااهمیت یکی از روش‌های پيشرفتة 
خرابی و هم‌چنین نعیین محدوده‌های خرابی. مورد 
استفاده قرار می‌گیرد. اگر ناحیةٌ خرابی به‌دست آمده از 
تحلیل سازه‌ها ( باشد. آن‌گاه در حالت کلی احتمال 
خرابی هر سازه پراساس تحلیل قابللت اعتماد در این 
)۱ 00۳ ]| - ,0 


/ 
در رابطة (۰۱ ,2 تابع چگالی احتمال مشترک 
متغیرهاء 1 ناحیهٌ خرابی معین شده به‌وسیلهٌ یک يا چند 
تابع حالت حدی و 3۶ میزان احتمال خرابی سازه 
می‌باشد. با تعریف تابع نشانگر («10 به‌گونه‌ای که به‌ازای 
نمونه‌های واقع در ناحيهء خرابی 10-21 و در سایر 
موارد 20 (*)1 باشد. می‌توان رابطه (۱) را برای همه 
متغیرها ب‌صورت رابطة (۲) بیان نمود [ ۲]. 


0, - | 10, )(۵ 3 


/ 
برای تولید و شبیه‌سازی نمونه‌هاء به‌جای استفاده 

از تابع توزیع احتمالی اصلی متغیرهاء («8,0 از یک تابع 
چگالی احتمال جدید (,» با عنوان تابع چگالی 
نموت گیری یاهمیت فاد بی‌شره تارمن یج 
شبیه‌سازی باید به‌صورت توابع وزنی برای تخمین تابع 
چگالی بااهمیت. به‌کار گرفته شوند. باتوجه به موارد 
ذکر شده. احتمال خرابی رابطهٌ (۲) به‌صورت رابطةٌ (۳) 


تصحیح می‌شود. 


/ 


3 
9 هرد ۳ < ,0 


در رابطة (۳) تابع (10۷ به‌عنوان یک تابع نشانگر 
می‌باشد. مقدار این تابع برای نقاطی که در محدوده 


سال بیست و نهم شماره یک ۱۳۹۲ 


۱۱ 


خرابی واقعند پرابر ! و در غیر این صورت برابر صفر 
درنظر گرفته می‌شود. تابع ()," نیز به‌عنوان تابع چگالی 
بااهمیت مناسب می‌باشد. یک تخمین نااریب از احتمال 
خرابی به‌صورت رابطة (4) بیان می‌شود. 


_ 1 02 
۳ 0 3 


در رابطةٌ (6) نمونة تصادفی :۷ بااستفاده از تابع 


توزیع بااهمیت (),۲ تولید شده و (] نیز تابع توزیع 
احتمال اصلی متغیرها است [15]. 


روش نمونه‌برداری ابرمکعب لا تینی ۲161 
((۳۲۲5) مصامرصجه امعم وو۲ متام 

روش نمونه‌برداری ابرمکعب لاتینی برای اولین بار 
توسط «مکی» در سال ۱۹۷۹ ارائه شد [16]. این روش 
براساس روش لایه‌لایه‌شده (۱6/00 1۵6۵۲60 ۲۳۵) به 
تقسیم محدوده تغییرات هر متغیر به‌صورت احتمال 
مساوی می‌پردازد. بدین ترتیب فضای متغیرها ماه 
نمونه گیری از محدوده هر یک از این ابرمکعب‌ها 
به‌عمل می‌آید و براساس قضية احتمال کل م۳ به‌دست 
می‌آید. در این روش تعدادی از مقادیر ممکن هر متغیر 
که به لایه‌هایی جزءبندی شده است. انتخاب می‌شود و 
یک مقدار از هر لایه به‌طور تصادفی به‌عنوان مقدار 
نماینده انتخاب می‌گردد. تابع حالت حدی ۷ با ک متغیر 
)0( ری ۲ 

برای تست اروذن اولین ترکیب از هر متغیر 
تصادفی کل یک نماینده به‌طور تصادفی انتخاب می‌شود. 
برای به‌خست آورون دومین ترکیب. یک نماینده به‌طور 
تصادفی از ۱۷-1 مقدار باقی ماندهٌ هر متغیر تصادفی > 
انتخاب می‌شود و این مراحل تابه‌دست آوردن ۲ 
ترکیب از متغیرهای تصادفی ادامه داده می‌شود و 
درنهایت [ا مقدار نماینده برای هر کدام از >1 متغیر 
تصادفی به‌دست می‌آید. در کل» ۷ ترکیب محتمل از 
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۱۱۲ 


این مقادیر می‌تواند وجود داشته باشد. همدف روش 
ابرمکعب لاتینی این است که با انتخاب آا ترکیب از 
این مقادیر نماینده» یکی و تنها یک مقدار از آا ترکیب 


تست و 


فرمول‌بندی روش رگراسیونی کوچک‌ترین مربعات 
دستگاه پردارهای پشتیبان (1,55*1 [10-11] 
تابع توزیع معین (۳)۶,۷ که در آن ٩۳,۷‏ > ۶« است. 
را درنظر بگیرید. درنتیجه می‌توان نقاط نمونه‌برداری را 
پراساس رابطة (1) به‌صورت نگاشت نقاط در فضای 18٩‏ 


به 1 تولید نمود. 


0( تن ی ولا وق ورب :1 
در رابطه (۲,۰6۱۱ بردارهای ورودی» ۷ مقدار 


,11,2 4 ابعاد بردار و 1 تعداد 


تابع حالت حدی؛» ا, 
مشاهدات ورودی‌ها و خروجی می‌باشد. می‌توان رابطءهٌ 
بین بردارهای ورودی ‏ و خروجی ۷ را به‌صورت 
(۶6- ۷ درنظر گرفت. اگر رابطهٌ خطی بین نقاط 
نمونه‌برداری. فرض شود آن‌گاه تابع رگراسیونی 
به‌صورت 9 +(۷*06 < (20 به‌دست می‌آید. در ایین 
تابع ۷۷ و 9 به‌ترتیب بردار وزن و عرض مدل به‌دست 
آمده می‌باشد. خحطای آموزشی تابع براساس کمترین 
ریسک به‌صورت رابطة (۷) به‌دست می‌آید. 


1 
2 ۵-۳۷۱ ,)ز«زده 


احز 
در رابطهٌ (0۷ 0  <‏ و تابع ‏ براساس روابطی که 
در ادامه ذکر می‌شود به‌دست می‌آید؛ اما در اغلب موارد 
ارتباط غیرخطی بین بردارهای ورودی و خروجی 
برقرار می‌باشد. درنتیجه تابع رگراسیونی براساس تابع 
غیرخطی ‏ به فضایی با ابعاد بزرگ‌تر تبدیل می‌شود و 
درنتیجه رگراسیون خطی به رگراسیونی با ابعاد بالاتر 
تبدیل می‌شود. درنتیجه اگر بین نقاط نمونه‌برداری 
ورودی و خروجی رابطهٌ غیرخطی (۸) فرض شود 
آن‌گاه روش 1,۹5۷ به‌وسیلةٌ تابع غیرخطی ۰ تابع 
رگراسیون خطی را به یک تابع غیرخطی در فضای 
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نمونهٌ اصلی تبدیل می‌کند. 


۲0,۷۷ ( < ۷۷ ۰)( + 0 (۸ 


بنابراین الگوریتم 1.85784 قادر به تعیین تابع 
رگراسیونی براساس یک مسئلهٌ بهینه‌سازی برمبنای 
برنامه‌ریزی غیرخطی (/|"|) طبق رابطة )٩(‏ و با قید 
تساو ۱ بش 


1 
9 7 ۷+ 2-۳ - ۵,")ز صرنصر 
1ح[ 


6 +9 +(06 ۵ ۷ :1۷ 
اس 1 


(۱۰ 


و با مشتق‌گیری از تابع لاگرانژ جمله‌ای قیود بهینه‌یابی 
نیز به‌صورت روابط (۱۲) به‌دست می‌آیند. 


1 
(ر- رو+0+( )رم - مزع (هیط م1 
[-ز 


)۱۱( 
1 
20 / رای ۳ 
1<1 00 
. 18 
3 2 2 0 اس 
6 2 ۳ تک 
,91 
< وب 0<س- 
1 1 1 :0 
20 - 04۵+( 0 اب 0-< 
0 
)۱۲( 


روابط (۱۲) را برای سادگی می‌توان به‌صورت 
ماتریسی به‌شکل رابطة (۱۳) نمایش داد. 


۲ 0 0 2-7 |۷7 | 0 
1 
0 0 0 1 9 5 0 ۳ 
0 0 ۷ -1 || 0 
7 1 1 0 2 ۷ 


(۱۶) به‌دست اند 
(۱( ۷ ی و 
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[(0 (یعت (رعه| 7 
ه 8 رعأ <ه 

قم 32:۰ ,2 < 2 

۳۹ 1 


باتوجه به این که در مسائل بهینه‌یابی محاسبهٌ 

ضرب داخلی (۵0,(۵/6 مشکل می‌باشد. بنابراین 

به‌منظور تسهیل در محاسبات می‌توان تابع کرنل 

(ر«,۲6 را در روابط (1۶) به‌جای ضرب داخلی نمود. 

(۱۵ رد۱۷6 ( ۵6 ۲ ((۵0) ع ی 8۵ 

)[<12....( 

برای تعیین پارامترهای تابع رگراسیونی انتخاپی (2 

و 0) می‌توان با حذف 6 و ۷ در رابطة (۱۶) و بااستفاده 
از وانظه (۱۵) رانطه (۱) را بفوست آوود: 


3 اه 


در رابطة (۱3) بافرض 1 +240 ۸ و باتوجه به 
معین بودن ماتریس فرض‌شده می‌توان پارامترهای تابع 
۱۷ ([0- 2-۸۱0 


بااستفاده از او تن رابطه از روابط (۱۲) و 


جایگزینی آن با رابطةٌ (۱۵) تابع رگراسیونی انتخابی 


به‌صورت رابطةٌ (۱۸) به‌دست می‌آید. 


۱۸( و و و 


الگوریتم‌های پیشنهادی براساس ترکیب 1,85۷ با 
روش‌های نمونه‌برداری بااهمیت و اپرمکعب لاتینی 
(۲۲9:]-۱۷ 1:5۷ [۲5]-15۱۷۲۷۲) 
ایده روش‌های پیشنهادی به‌دست وف مدلی کارا و 
مناسب برای توصیف مجموعه داده‌های آموزشتی و 


تخمین دقت مناسب برای داده‌های آژشایشین می‌باشد. 


سال بیست و نهم. شمارةٌ یک ۱۳۹۲ 


۱۳ 


برای ایجاد چنین مدلی باید براساس میزان‌سازی بین 
مدل آموزش‌دیده و مقادیر تابع حالت حدی واقعی, تابع 
پاسخ مناسبی را به‌دست آورد این در حالی است که 
مدل‌های دیگر باتوجه به عدم میزان‌سازی ممکن است 
توابع نامناسبی را به‌دست آورند. باتوجه به این که حل 
انتگرال مونت‌کارلوی معمولی رابطة (۲) در سازه‌همای 
خرپایی و قاب‌های فولادی که دارای توابع حالت حدی 
با متغیرهای زیادی هستند. نیاز به زمان بسیار زیادی 
خواهد داشت. بنابراین برای تسهیل در انتگرال گیری. از 
روش‌های نمونه‌برداری بااهمیت و ابرمکعب لاتینی 
استفاده شده است. با توجه به این که واریانس مقدار 1۲ 
باون6 (واریانس متغیرها) رابطه مستقیم و با تعداد 
نمونه‌های تولیدشده رابطةٌ عکس دارد درنتیجه انحراف 
معیار مقدار ۲ و حاصل تخمین مونت‌کارلو با نسبت 
۱۷ کاهش می‌یابد و این حکایت از همگرایی کند روش 
مونت‌کارلو دارد. همگرایی کند روش مونت کارلو با 
روش‌های 15۷ و 1۳15 که برمبنای کاهش واربانس 
هستند جبران خواهد شد. باتوجه به این که روش 
نمونه‌برداری بااهمیت برای افزایش تعداد نقاط 
نمونه‌گیری در ناحيةٌ خرابی. از یک تابع چگالی احتمال 
که تأکید آن برروی ناحية خرابی است. استفاده 
می‌نماید. روش مناسبی برای نمونه‌برداری داده‌ها 
می‌باشد. برای انجام نمونه‌گیری و پخش نمونه‌ها در 
قسمت‌های مختلف خرابی به شناسایی قسمت‌هایی که 
چگالی احتمال در آنها بزرگ است. اقدام می‌شود. روش 
نمونه‌برداری ابرمکعب لاتینی نیز برای تحلیل قابلیت 
اعتماد سازه‌های پیچیده‌ای که زمان محاسبهة فقط یک 
نمونة شبیه‌سازی آن ممکن است بسیار زیاد باشد 
به‌منظور کاهش تعداد نمونه‌های مورد نیاز مناسب است. 
باتوجه به مطالب ارائه‌شده براساس روش‌های پيشرفتة 
مونت کارلو تعداد تحلیل سازه نسبت به روش مونت 
کارلوی معمولی به‌طور قابل‌توجهی کاهش می‌یابد اما 
هم‌چنان تحلیل المان محدود امری ضروری می‌باشد. 
همان‌طور که در قسمت‌های قبل ذکر شد. روش‌هایی از 
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۱۹ 


فا که یی تالا سل له ستاوی تس 
تجربی و براساس توپولوژی مناسب تنها می‌توانند 
خحطای آموزش را کاهش دهند و تعميم‌دهنده کل داده‌ها 
نمی‌باشند و درنتیجه روش تحلیلی مناسبی برای 
پیش‌بینی رفتار سازه نمی‌باشند» این در حالی استت. که 
ساوبای :و دون تیان به تویرلرزی با بجاد عدلی تسب 
و تعمیم آن به کل داده هاء روش تحلیلی مناسبی برای 
پیش‌بینی رفتار سازه ارائه می‌ دهد و علاوه بر به‌دست 
آوردن خطای مناسب» زمان محاسبات را نیز به‌مقدار 
قابل توجهی کاهش خواهد داد. 

در اين مقاله الگوریتمی براساس ترکیب روش 
بااهمیت و ابرمکعب لا تینی ارائه شده است. اصول کلی 
داده‌ها که براساس روش‌های نمونه‌برداری بااهمیت و 
ابرمکعب لاتینی ایجاد شده‌اند. به دو گروه داده‌های 
آموزشی و آزمایشی تقسیم می‌شوند و پارامترهای مدل 
بااستفاده از آموزش مجموعه داده‌های آموزشی (آموزش 
مدل 1/58۷۱۷) و براساس توابع مختلف کرنل. به‌دست 
جدیدی محاسبه می‌شود و درنهایت مقادیر تابع 
رگراسیونی به‌ازای داده‌های نمونه‌برداری شده محاسبه 
می‌ شود و میزان احتمال خرابی تنخمین زده می‌شود. 
الگوریتم‌های ترکیبی براساس اصولی که هر کدام از 
روش‌ها دارند. برتری‌های خاصی نسبت به تک‌تک آنها 
دارد. یکی از این مزیت‌ها. کاهش تعداد نمونه‌ها توسط 
روش‌های نمونه‌برداری پیشنهادی و کاهش تعداد 
تحلیل‌های المان محدود برای تعیین احتمال خرابی 
است که در صورت استفاده از مدل گر اتسیو نی 
روش‌های نمونه‌برداری برای تعیین نمونه‌ها و مدل 
۲ 1[ برای تقریب محدوده خرابی به منظور تسهیل 


ذر انتگرال گیری عددی به‌کار می‌رود. در ادامه به بررسی 
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تحلیل قابلیت اعنماد سازه‌ها براساس روش ت رکیب ی کوجک‌انرین... 


تحلیل قابلیت اعتماد سه سازه با توابع حاللت حدی و 
تعداد متغیرهای مختلف پرداخته و کارایی بالای 
الگوریتم‌های پیشنهادی نسبت به دیگر روش‌های نوین 
نشان داده می‌شود. تعداد داده‌های آزمایش برابر ۲۵ و 
تعداد تحلیل‌ها برابر ۵۰ درنظر گرفته شده است. نتایج 
نشان‌دهندة کارایی بالای روش‌های پیشنهادی نسبت به 
دیگر روش‌ها می‌باشد. به‌منظور تشکیل مدل کارا برای 
هر سازه. تحلیل حساسیت بر روی مجموعه داده‌های 
آموزشی انجام شده و مناسب‌ترین تعداد داده برای 


افق دی بازفیدکه املع ابیگ: 


مثال ا: قاب پرتال فولادی با ۳ دهائه و 4 طفه. 
تحلیل این قاب برای اولین بار توسط آقایان بوچر و 
بورگاند ارائه شد [15]. شمای این قاب در شکل (۱) 


نشان داده شده است. 


شکل ۱ شمای قاب پرتال حطی ۳ دهانةٌ ۵ طبقه 


این قاب پرتال شامل ۲۱ متغیر تصادفی می‌باشد. 
این متغیرها شامل مدول یانگ ممان اینرسی؛ سطح 
مقطع اعضا با توزیع نرمال و هم‌چنین بارهای جانبی با 
تابع توزیع رایلی می‌باشد. خصوصیات و پارامترهای 
آماری این متغیرها براساس جداول (۱) و (۲) می‌باشند. 
اگر جابه‌جایی مجاز این طبقه برابر با ۱۰ سانتی‌متر 
باشدء آن‌گاه معیار خرابی براساس تحلیل المان محدود 
به‌صورت تابع حالت حدی رابطة (۱۹) می‌باشد. 

0060-10-1, ۱( 
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در واه (۵ یلا اب بای اففتی وافنین 
برحسب تابعی از متغیرهای تصادفی می‌باشد. نتایج 
به‌دست آمده براساس الگوریتم‌های پیشنهادی هم‌چنین 
نتایج روش‌های مراجع دیگر در جدول (۳) ارائه شده 
رمفت )یال شش اس دس رشان روش 
شبیه‌سازی مونت کارلوی معمولی و با ٩۰۰۰۰۰‏ نمونه 
بعة از *۱ باز تخلیل پیانی در زمان ۹۵۱۲ کائيه برابر نا 
۹ به‌دست آمده است). نتایج ب#دستت: اسلا ون 
لول (۳) نش ان‌دهتد۶ کاراین الا وشن هاق 
پیشنهادی از لحاظ دقت محاسبات و کاهش هزينة 


تخاسیات بت به دیگر روش‌هتای تعلیتل آیتن,قتات 


جدول ۲ پارامترهای آماری قاب ۳ دهانهة ۵ طبقه 


واحد انحراف معیار 


لک 9 


604 
41000 
1000 

13-۳207 1" 
1۱-1 
3-14 
018 
3-14 
+352 
1319 
2-2-33 
6۰157 
3 
10692 
1-314 
109 
162 
13.900 
1۱-۰304 


۷-۷-۱-0 


۳ 


تفت : 


میانگین 
30 
20 
16 
5400( 
19100 
10۳0۳۰7 
133206 
2-15 
2-14 
142.92 
19-5۰-3711 
2-69 
2-2-14 
2-1 
16( 
22-1093 
62۳-1534 
204 
2368 
1(2978 
905( 


سطح مقطع 


جدول ۱ خصوصیات اعضای قاب ۲ دهانة ۵ طبقه 


۳۳ 
۸6 
۳ 
۸8 
۸1 
۸2 
۸3 
24 


تلع توزیع 
۹2/12۲ 
۹2/12 
ات۸ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 


۱۹ 


۱۹ 


ممان اینرسی 
15 
16 
17 
18 
11 
12 
13 
1 


متغیر تصادفی 
۳1 
۳2 
۳3 
1۳1 
۳2 
11 
12 
13 
14 
15 
16 
17 


جدول ۳ نتایج روش‌های پیشنهادی و مرجع در مثال ۱ 
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مدول یانگ 


1۳1 
1۳1 
1۳1 
1۳1 
۳2 
۳2 
۳2 
۳2 


اعضا 
181 
۳2 
۳3 
۳94 
اتف 
42 
98 
4 
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۱۹3۹ تحلیل قابلیت اعتماد سازه‌ها براساس روش ترکیب ی کوچک‌ترین.. 
خحطای نسبی (/) | زمان (ثانیه) | احتمال خرابی تعداد داده آموزشی | تعداد نمونه روش 
0342 142 090029 ۹00000 1,5۷5 
092361 2467 0000291 200 /- 1/551۷ 
1.02397 1۹11 0000292 20 230 ۱92 
2252 02 1 0000291 : [8] ۳0۳0/۲ 
63 3 ۰ 000022 ِ [8] 112101۳760 
78 ِ 02000312 روش سطح پاسخ کلاسیک [8] 


مثال ۲: قاب پرتال فولادی با ۳ دهانه و ۱۲ طیقه. 
قاب مفروض شامل یک قاب پرتال الاستیک خطی با ۳ 
دهانه و ۱۲ طبقه می‌باشد [7,13]. شکل (۲) این قاب را 
به‌همراه تیپ‌بندی تیر و ستون‌های آن نشان می‌دهد 
(ارتفاع هر طبقه 4 متر می‌باشد). سطح مقطع تیر و 
ستون‌ها و هم‌چنین بار افقی وارد به قاب به‌عنوان 
متغیرهای تصادفی مستقل درنظر گرفته شده‌اند. مدول 
یانگ به‌عنوان پارامتر قطعی برابر با "۲۲۱۰ می‌باشد. 
ممان اینرسی اعضای این قاب نیز به‌عنوان متغیری 
وابسته به سطح مقطع‌هاء براساس رابطه (۲۰) بیان 
من 3 
(۲۰( توا[ 
با فرض این که ماکزیمم جابه‌جایی افقی گره ۸ 
برابر با ۰/۰۹7 متر باشد. آن‌گاه تابع حالت حدی 
براساس پارامترهای تصادفی فرض شده به‌صورت 
رابطة (۲۱) درنظر گرفته می‌شود. 


(۲۱( (۳ وی ,۸۵ و۸۵ ,یط ,)را -0.096 2 (000 


پارامترهای آماری نمونه‌های تصادفی اصلی و 
ضرایب ممان اینرسی این قاب نیز مطابق جدول (۶) 
می‌باشد. نتایج به‌دست آمده از روش‌های پیشنهادی 
براساس تعداد نمونه‌های آموزشی که براساس تحلیل 
تخباسیت ب‌دست امد و ادله به‌همراه نتایج روش‌های 
پیشنهادی مرجع اصلی به‌صورت جدول (۵) می‌باشد. 
نتایج به‌دست آمده نشان می‌دهند که روش‌های 


ی پیشنهادی هزینة زمانی را به نحو موثری کاهش می‌دهند 
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و برآورد احتمال خرابی را با دقت مناسب ارائه می‌دهند 
(مقدار احتمال خرابی دقیق با ۱۰۰۰۰۰ نمونه 
مونت‌کارلو و با ۱۰ بار تحلیل پیاپی در زمان ۵۲۶۱ ثانیه 
برابر با ۰/۱۷۵۱ به‌دست آمده است. این میزان برابر 


مقدار دقیق ارائه‌شده در مر جع ۷ می‌باشد). 


شکل ۲ شمای قاب ۳ دهانةٌ ۱۲ طبقه 


جدول ۶ پارامترهای آماری قاب ۳ دهانه ۱۲ 


سال بیست و نهم. شمارهٌ یک ۱۳۹۲ 


محمد مهدی مجاهد- بهرو احمدی ندوشن 


اون 5 انحراف 7 ۲ 
1 توزیع واحد »۳ میانگین | متغیر 
3 0و1[ ر 0025 025 ۸۱ 
3 0و0[ 9 0016 016 ب۸ 
3 و1 0036 036 وه 
026617 موصعم [ 9 002 02 بر 
02 ۱ ر 0015 015 و 

2 30 75 ۹ ِِ ِ 


مغال ۳ خرپای ۱ عضصوی با تایع حالست حدی 
جابه‌حایی. مثال آخر یک خریای ۱۵ عضوی است که 
از مرجم ۳ برگرفته شده است. هندسة این خرپا در 
شکل (۳) نشان داده شده است. سطح مقطع همه ۱۵ 
عضو و ۳ بار متمرکز که بر گره‌های ۲ ۵ و ۱ وارد شده 
می‌باشند. مدول یانگ نیز به‌عنوان پارامتر قطعی برابر با 
۸ ۰ فرض می‌شود. براساس تحلیل الاستیک 
خطی, میزان جابه‌جایی مجاز گرة میانی (گره ۵) برابر 


با ۷/۵ سانتی‌متر به‌دست آمده است. پارامترهای آماری 


۱۱۷ 


متغیرهاء در جدول (1) ارائه شده است. هم‌چنین نتایج 
به‌است امه از روتن‌شضای, پشتهاذی بر اساس. نداد 
نمونه‌های آموزشی که براساس تحلیل حساسیت 
به‌دست آمده‌اند به‌همراه نتایج روش پیشنهادی مرجع 
اصلی به‌صورت جدول (۷ می‌باشد. میزان احتمال 
خرابی دقیق در مرجع ۲ برابر ۰/۰۱۶ به‌دست آمده 


است. 


ات 90 4.6 سس 10 6 4ص (0] 4.6+( 100 4.6 


شکل ۳ هندسه خرپای ۱۵ عضوی 


جدول ۵ نتایج روش‌های پیشنهادی و مرجع در مثال ۲ 


خطای نسبی | زمان (ثانیه) | احتمال خرابی 
0.065 1109 00515«( 
0.0931 57 2517( 
0.65 5.15 20515«( 
6 1917 00115 0 
2.346 5907 0.0117 
2-417 ۳ 00405 


امصصم 9 
امصصمآ تن 


امصمتع][ الک 
امصصمتع ] لک 


امصمتع0] الک 


تعداد داده | تعداد نمونه روش 
100 100000 ۱۷-56 1,85۷ 
00 1000 /7-]۱۷ 15۷ 
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انحراف معیار | میانگین متغیر 
0316 1032 ۸1-6 
0393 6045 ۸7-5 
۰09975 0 3 + 1 1 
#09975( 382 2 1 
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جدول ۷ نتایج روش‌های پیشنهادی در مثال ۳ 


سال پیست و نهم شمارة یک. ۱۳۹۹ 


نشریة مهندسی عمران فردوسی 


۱۸ 


تحلیل قابلیت اعنماد سازه‌ها براساس روش ت رکیب ی کوجک‌انرین... 


خطای نسبی (/) | زمان (انیه) | احتمال خرابی | تعداد داده های آموزشی 
5- 1936 2( 50 
4- 2.92 95( (0 50 
94- 375 2.139( 00 
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تحلیل حساسیت برروی داده‌ها ی آموزشی و تعیین 
میبزان خطای داده‌های تست در مشال ۲. تحلیل 
حساسیت بر روی مدل. براساس تعداد مختلفی از 
داده‌های آموزشی و با ثابت درنظر گرفتن تعداد داده‌های 
تست (تعداد داده‌های تست برابر با ۲۵ درنظر گرفته 

شده است) انجام می‌شود. براساس این تحلیل 
حساسیت. محدوده مناسبی از داده‌های آموزشی به‌دست 
می‌آید. همان‌طور که در شکل (4) قابل مشاهده است؛ 
در هر دو روش پیشنهادی ۲۲5,]-95۷۱] و 85۷۷۲,] 
میزان جذر میانگین خطای مربعات (130718۳), 
ناشی از مقادیر تابع حالت حدی مدل ایجاد شده برای 
داده‌های تست باافزايش تعداد داده‌های آموزشی. 


تحلیل حساسیت برروی داده‌ها ی آموزشی و تعیین 
میزان احتمال خرابی در مثال ۲. در این قسمت. تحلیل 
حساسیت برای یافتن بهترین مجموعه از داده‌های 
آموزشی. براساس میزان احتمال خرابی به‌دست آمده 
انجام شده است. این تحلیل بدین صورت است که در 
هرکدام از روش‌های پیشنهادی, با درنظر گرفتن تعداد 
مختلفی از داده‌های آموزشی, میزان احتمال حرایی سازه 
براساس ۱۰۰ بار تحلیل به‌دست ی رش نتایج به‌دست 
آمده در جداول (۸) و )٩(‏ نشان داده شده است. شکل 
(۵) نیز مقادیر احتمال خرابی را برای تعداد داده‌های 
آموزشی مختلف برای هر دو روش نشان می‌دهد. مقدار 
ضریب تغییرات کمتر از ۰/۱ و در صد خحطای کمتر از 
یسک درصد ببه‌عضوان مقادیر مناسب درنظر گرفته 
می‌شوند. نتایج ب4‌قستتا آبله نشان می‌دهد که با افزایش 


تعداد داده‌های آموزشی. میزان ضریب تغییرات و درصد 


خطای نسبی مقادیر احتمال خرابی کاهش می‌یابد. 
باتوجه به این که مقادیر ضریب تغییرات و درصد خطای 
احتمال خراپی به‌دست آمده براساس تعداد ۵۰ و ۷۰ 
داده بهترتیب برای روش‌های ,9۷۱ص و 
۱۷-۲ 9۱۷ در بازه مناسبی قرار دارد. درنتیجه 
می‌توان براساس تحلیل حساسیت داده‌های آموزشی, 
این تعداد را به‌عنوان تعداد داده‌های آموزشی مناسب 


درنظر گرفت. 


جذر مبانگین خطای مربعات 


د فا 
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تعداد داده های آموزشی در روش ۱ 


شکل ۶ نمودارهای جذر میانگین خطای مربعات براساس تعداد 
داده‌های آموزشی مختلف 


جدول ۸ تحلیل حساسیت برروی داده‌های آموزشی و تعیین احتمال خرابی در روش 1115-/1.55۷ 


نشرية مهندسی عمران فردوسی 


سال بیست و نهم. شمارهٌ یک ۱۳۹۲ 


محمد مهدی مجاهد- بهرو احملی ندوشن 
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جدول ٩‏ تحلیل حساسیت برروی داده‌های آموزشی و تعیین احتمال خرابی در روش 153-//1.85۷ 
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تعداد داده های آموزشی در روش /(1,88۷۰13 


شکل ۵ نمودارهای میزان احتمال خرابی براساس تعداد مختلفی از داده‌های آموزشی 


جع پندي 


سال بیست و نهی شمار؛ُ یک» ۱۳۹۲ 


احتمال خرابی براساس ردش ,13512-193 


۱۹ 


نشریة مهندسی عمران فردوسی 


و ]151-]1858۷ پرداخته شد. براساس این روش‌ها و 
بااستفاده از داده‌های آموزشیء مدل مربوط ایجاد شد و 
درنهایت براساس روش‌های پیشنهادی. میزان احتمال 
خرابی تخمین زده شد. با ارائة ۲ مثال از قاب و یک 


مثال خرپایی. نشان داده شد که روش‌های پیشنهادی 


تحلیل قابلیت اعنماد سازه‌ها براساس روش ترکیبی کوچک‌نرین... 


قابلیت اعتماد نیاز به زمان زیادی جهت تحلیل المان 
محدود خواهند داشت. هزینهة زمانی صرف‌شده از طریق 
روش‌های پیشنهادی به‌طور قابل توجهی از روش‌های 
نکر کم می‌باشد. در ضمن. دقت محاسبات نبز در 


حد معقول و حتی بهتر از روش‌های دیگر خواهد بود. 


جهت تقریب تابع رگراسیونی. برای توابعی بسیار مفید 
هستند که دارای فرم ریاضی نیستند و جهت تحلیل 


ترا 
2000(۰) و1-6 .00 ر01)ک0ظ .۱۷66۲۵ و قمتنطامنای ۵۴ وطتلاطافتامط ریک رفصتلامن روم باه 


م۳ عصمتامص۲ ماهاگ )نصا اتمتامرا مصتامند۸۵0 1۵۲ قلمطام]۱ ۵۶ صمت)مصتصصید۲ صه ویظ.ژ ومه) ۲۳ 
(2004) ,271-293 00۰ ,26 ۷۵۱۰ ,029 ایک ۲3۵۵۲۲ عصتصنجما آهمااوتاماگ 0۴ ۷۱۵۷/۵۵ عطا 
مناد 10۲ ولمطاع]۱۷ مایمن مم۵ [72مص1 ۲ رب رمه‌تعصاه 0صه ری رقصصهط] ربا راته00ظ وب رع7۳0208 
(2010) ,183-190 00۰ ,32 ,۷۵۱ ,50/69 آم امک و انازطادنای؟ 

ولافک صعتوعض 20۷۵۲۱ 2۲ مفطمتتوظ روتاامفتامف؟ اجستمند ظ ومممه. معممم‌ومط ویک ۲۹۵26017162 
1962(۰) مصعطمص ۱۷۲ 67 ۲۱۵۲۲ م77 20ظ هل 1۵0118[ 

۵ 01امرحصا بو وتورافمم وتااممتاف؟ امستمناو 72.0 ر۵ ۷۲ 200 ,2 ون ولا بتات) و.ل ,16182 
.(2005) ,25-48 .0 و27 ۷۵۱۰ 6۵/69 ماگ ۱۱۵۱۲۵۲ هبو م۸1 عصونا مصمتامعصا۲ 

مان ۲۷۷1 ۱۱۵۱۷۵۱1 منوا مج وموکو مقصممروم ۵۶ ممواهم۵ن .ریخ مه .۲۱ رقععصمت) 
2004(۰) ,49-67 .00 ,26 ,۷۵۱ نواعم میتی رقلو همم الااطامتای؟ ماگ بم۶ م0مطاه]۱۷ 

0 0۳۵۵۵5 کل امصه متا لمتمناد ‏ به ‏ فصطاتدموله ملامصمت ۳۳۱۷/۵۳۵۵ .لگ رصق 
.(2007) ,1524-1533 .00 ,85 ,۷۵۱ بکه50۳16 

,6 ایک رقصهمامام۳ هناگ لمیتمناوه 6۵ ۱۷۲۵0۵0 مصتع‌تنک ۵۶ ممتاهءلآ۸۵ ریگ 1521۷۳222 
(2005) ,133-151 00۰ و27 ,۲۷۵۱ 

1998(۰) ۷۲۵۲ و۱۱ ر.عط0گ 24 ۷۷۱۱۵ رهز ۲6۵0۲۲ مصتصریهمنا اممتاعتاه۹" ۲۷۰ ۲۷20۵۲1 

0 ۵۵۲۵۹۹100 10۲ مصتطمه۱ ۷۵۵۲۵۲ ا۲ممرمرتا٩‏ ومتهنووی اقععا ۵۶ مملاهعنآممرظ" ر.ظ رصهطه‌عصهنان رت رل70107 
.(2009) ,160-166 .00 22 ,۷۵۱ رکع۵۲۸۵۲0۳۵ ول ۱۱۱۵5۵ ولو ومصم الاامامنای؟۲ 

0 ,"۸217۹15 واتااطمتای؟ مناد عم ومتوومتوم؟ زما۵۵ ۷ ]۵0۲صتا9 و2 وتان ۸0۵ .1.9 بان 
.(2006) ,1295-1303 .00 و(27)10 ۷۵1۰ وکن۷66۳۵( ۵060 ۷۲۵۲۱۱6۱۱۵۵66( 


هام0 ما مملامعنآمروج کل 240 0مطامصطر ممهتاه مقممموع۲ کصماماگه صه و2 بتالم ,۷۷ ,2020 . 


,34-42 .00 و67 ,۷۵۱ متواعع۲ ۵80 کادرا2۳۵ 1۱ ۱6۳666 مانا رصمتاحعتصتانن 0عفوها- ماع 20 
(2013) 
65 ۲۵0۱21 عصنوعده صمتاعصا۲ ممصهمصوم)۳۵ )امتام‌حیا ت۶۵ فلورافمم وتلاطامتاف؟ لقتنتمتتطو رز رعصهرز 


(2006) 3255-3291 .00 ,43 ,۷۵۱ رکع نت60 وه کلزآوه ]۵ ماو امص۲۱۵۵۵0 ۵۱۷۵۲۱ ممتامصتا۲ 


نشریة مهندسی عمران فردوسی 
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11۰ 


13۰ 


سال بیست و نهم. شماره یک ۱۳۹۹ 


۲ ۲0۳9۳۶ رز 5۷۱۷۲ 0۶ 60وده ۱۲۵/۵0 عصتام‌صحه عومصماعممرصا ر.گ رصعطن 20 ۷۷.۲ وطال ویک ر1218" .14 
2007(۰) ,62-5 .00 ,(24)6 ,۷۵۱ ما75 5016۳61۶65 

امستمنی عم طمهمتووخ ممعسته معممموم؟ اصملم ۳۲ صرح بعه مق ولا مصاعتتا0ظ .0 متعطهناظ ,15 
(1990) ,57-66 .00 و(701 ۷۵۱۰ امه تناو رعصمامامظ واتلامامتای؟] 

توص که مماه ۲۷ عصتامعام ۶۵۲ ولمطاع۱۷ عع۲ظ۲ ۵۶ جموتنعم‌ههه خر و.لیظ؟ م۲601 رای و۱۵۵ .16 
,239-45 .00 و(42)1 ۷۵۱۰ 16611۳00۳1665 ,"00۵06 میامن 2 حصم کبانان ۵۶ فلووآهمه مط صا ومامام2۳1 ۲۷ 
2000(۰) 


سال بیست و نهم. شمارةٌ یک ۱۳۹۲ تشن هنم :عمران فردوسی 


۳9۲ تحلیل قابلیت اعتماد سازه‌ها براساس روش ت رکیب ی کوچک‌نرین... 


نشریة مهندسی عمران فردوسی سال بیست و نهم. شماره یک ۱۳۹۱ 


